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ABSTRACT 



The formal theory of breakup reactions is reviewed. The direct breakup mecha- 
nism which is formulated within the framework of the post form distorted wave 
Born approximation is discussed in detail. In this theory, which requires the 
information about only the ground state wave function of the projectile, the 
fragment-target interactions are included to all orders while fragment-fragment 
interaction is treated only in the first order. The general applicability of this 
theory to describe the breakup of halo nuclei is demonstrated by comparing the 
calciilations with data for total, as well as energy and angle integrated cross 
sections and momentum distributions of fragments in reactions induced by a 
number of halo nuclei. We investigate the role played by the pure Coulomb, 
pure nuclear and the Coulomb-nuclear interference terms. Postacceleration ef- 
fects in the Coulomb breakup of neutron halo nuclei is also studied. 
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1 Introduction 



Recently advances made in accelerator and related technologies have provided 
us the opportunity to produce and work with nuclei having very short half-lives 
and very small one- or two-nucleon separation energies [1-14]. These nuclei 
lie very close to drip lines (the limit of neutron or proton binding). Nuclei 
at extremes of binding can exhibit behaviours which are quite different from 
those of the stable isotopes. We still lack a fully microscopic understanding of 
the stability of these unique many body systems. These nuclei are important 
also from the nuclear astrophysics point of view. The rapid neutron capture 
(the r-process) together with the slow neutron capture (the s-process), which 
are dominant mechanisms for the nucleosynthesis of heavy elements above iron 
pass mostly through the neutron rich region. The properties of these nuclei are 
important inputs to theoretical calculations on stellar burning which otherwise 
are often forced to rely on global assumptions about nuclear masses, decays and 
level structures extracted from stable nuclei. 

The first round of measurements involving neutron rich nuclei [151 1161 117L 
I18L I19L on] have confirmed the existence of a novel structure in some of them in 
which a low density tail of loosely bound neutrons extends too far out in the 
coordinate space as compared to a stable core (also known as the neutron halo"'^). 
The quantum mechanical tunneling of very loosely bound valence neutrons leads 
to the formation of such a structure. The existence of neutron halo has been 
confirmed in "Be ^''B [23 HI, ^^C [53] (one-neutron halo), and ^He, "Li 
Uni 1101123, "Be EHllini, and ^^B ^ (two-neutron halo). Recently, some 
proton halo nuclei have also been identified (*^B [301 EH El, ^^Ne [SO], ^OMg 

[Oa, and 26,27,28p 

Halo nuclei, in most cases, have only one bound state (the ground state) 
and a broad featureless continuum. Thus, methods of conventional nuclear 
structure studies, namely, measurements of energies and spin-parities of excited 
states are not applicable in these cases. However, due to their small binding 
energies, they can be easily excited above their particle emission thresholds. 
Hence their breakup reactions in the Coulomb and nuclear fields of the target 
nuclei could be useful tools to investigate their structures. 

To be able to extract reliable structure information of halo nuclei from the 
breakup data, it is quite desirable to have a theory of these reactions which (1) 
is fully quantum mechanical; (2) treats the Coulomb and nuclear breakups as 
well as their interference terms consistently on an equal footing; (3) includes the 
recoil of the core within the halo nucleus, and the finite range of the core- halo 
interaction, and (4) involves least adjustable parameters. 

We discuss the formal theory for halo breakup reactions in section 2. We 
specially mention a theory formulated within the framework of the post form 
distorted wave Born approximation (DWBA) where both Coulomb and nuclear 
breakups can be treated consistently on an equal footing. The full ground state 

^The term 'neutron halo' was introduced by Burhop et al. [211 . in the context (now called 
'neutron skin') of the bulk of the neutron density extending further out in space than the 
proton density. 
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wave function of the projectile enters as an input into this theory. Thus, infor- 
mation about the halo structure can be extracted by comparing the calculations 
with the available data. This is shown in section 3. We also discuss the role 
of nuclear breakup and of Coulomb-nuclear interference (CNI) terms in this 
section. The post form DWBA theory is uniquely suited to study the postac- 
celeration effects in the halo breakup reactions. This is an higher order effect 
which is studied in section 4. We give conclusions and the future outlook in 
section 5. 

2 Formal Theory of Breakup Reactions 






Figure 1: The direct breakup mechanism in a participant - spectator model. 

The basic mechanism of the breakup reactions can be described in a sim- 
ple participant-spectator model in which the projectile a which is supposed to 
consist of two substructures, say, b and c, interacts with a target t. It might 
so happen that one of the substructures b (the spectator) misses the target and 
keeps moving in its original direction while the substructure c (the participant) 
interacts with it. In such a situation the velocity of fragment b (vb) can be 
written as 

Vb = Va+Vp, (1) 

where is the velocity of projectile a, and vp is the velocity associated with 
the Fermi motion of b inside the projectile before its breakup (see Fig. 1). 
Therefore 

= vl+v% + 2v„.VF . (2) 

Asvp << w^, we can neglect Vp in Eq. (2) and get the maximum and minimum 
values of (corresponding to + and — signs, respectively) to be given by 

vl^vl±2vaVF. (3) 
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Hence, maximum and minimum possible energies of b in the final channel are 
given by 



f 2 mi) / 1 2 I \ / A\ 

H = 7:mbV^ = (o™a^a ± maVaVF), (4) 

2 nia 2 

where is the mass of fragment 6, and rua is the mass of the projectile a. 
Defining pa = maVa and pp — m\,vp, we have 

E,= '^Ea±P-^. (5) 
mo ma 

Thus, in the energy spectrum of particle b one would expect a peak at ^Ea, 
with a width (2 ^^-^ ) depending on its Fermi momentum inside the projectile. 
This simple picture, which was proposed by Serber |36) way back in 1947, is 
realized remarkably well in majority of the breakup data. 

Now, unless some kind of exclusive measurement is made in which both 
fragments b and c are detected, the fragment c can interact with the target 
nucleus, t, in all possible ways. We distinguish between two classes of breakup 
processes: elastic breakup, where the target nucleus remains in the ground 
state (the interaction c — t is elastic) and inelastic breakup where c can interact 
in all possible ways with the target nucleus (inelastic excitation of the target, 
compound nucleus formation, transfer processes etc.). 

There is another picture of breakup reactions in which the projectile is ex- 
cited as a whole to a state in the continuum which subsequently decays into 
fragments b and c as it leaves the interaction zone. This process is called the 
sequential breakup (or resonant breakup). Both Coulomb and nuclear interac- 
tions between the projectile and target can cause the inelastic excitation of the 
former. 

In the next sub-section, we shall describe the transition amplitude and its 
various representations and introduce the distorted wave Born approximation. 



2.1 The Transition matrix and distorted wave Born ap- 
proximation (DWBA) 

Let us consider the reaction a + t ^ & -|-c-|-t, in a three body model, where 
the projectile a, incident with momentum q^, breaks up into fragments b and 
c with momenta qf, and qc, respectively in the Coulomb and nuclear fields of a 
target t. 

The Hamiltonian of the system is written as 

H = n + T, + Tt + Vb, + Vu + Vcu (6) 

where Ti is the kinetic energy of particle i and Vij is the two-body interaction 
between i and j'; their separation will be denoted by in the following. 

To find the interaction in the initial and final channels, we note that the 
asymptotic Hamiltonians in the initial (prior) and final (post) channels are 

i/, = Tb + + (7) 
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and 

Hf=n+T, + Tu (8) 
respectively. Hence the initial (prior) channel interaction is 

V, = H-H,^ Vu + V,t (9) 
and the final (post) channel interaction is 

Vf = H - Hf ^Vbc + Vu + Vcf (10) 
There are two exact T - matrices [57] 

and 

^ ^^(-)|^^^ ^ Kt|e''>"""''^a(r..)), (12) 

which are the starting points for a discussion on the theory of breakup processes. 
The ground state wave function of the projectile, (j)a{vhc), satisfies 

(Tb + + Vbc)^a{rbc) = -CaMnc), (13) 

where Ca is the separation energy between fragments b and c in the ground state 
of the projectile. is the exact three-body scattering wave function with 

outgoing wave boundary condition (denoted by (+) sign in the superscript), 

and is the exact three-body scattering wave function with ingoing wave 

boundary condition (denoted by (-) sign in the superscript). They are the exact 
eigenfunctions of the three-body Hamiltonian [Eq. (6)]. Thus they satisfy 

= (14) 

and 

H^Y^^E^\r\ (15) 

where E is the total energy of the system. We now use the Gell-Mann - Gold- 
berger two potential formula [211 to rewrite Eqs. (11) and (12) as 



and 



= {x[-Av,t)x[-Ant)\v,. + + - - c/ct|*^^) (i6) 



In Eqs. (16) and (17), wave functions x[c \^ct)x[b \^bt) and x^t^\^at) satisfy 
the Schrodinger equations 

[Tr.,+Tr,, + Uu + t/,^]x(:Hrct)x(:)(rbt) = Ex[:\v,t)xt\vu) (18) 



6 



and 

[Tr^. + Uat]xt\'^at) = {E + ea)4tH^at), (19) 

respectively. In Eqs. (18) and (19), Uu are auxiliary potentials acting between 
particle 'i' and the target. 

Assuming Vm = Ubt and Vd = Uct in Eq. (16), we have P^HO) 

= (x(:'(rc*)xl:Hr,,)|F,.|*i+)). (20) 

Let us now introduce the distorted wave Born Approximation (DWBA) [JT] 
for the exact wave functions and in Eqs. (20) and (17). 

If one assumes that the inelastic excitations of the projectile are small, then 
the wave function can be approximated by 

*!^^«x[:^(raO0a(r.e). (21) 
The post form DWBA T - matrix is then 

tW^vosA^DwBA) = {x[-\vu)xtHvct)\V,Ax[':\vaMa{v,c)). (22) 

If, on the other hand, one assumes that the final state interaction between 
the breakup fragments {h and c) is not important, i.e. Vbc is weak in the final 
channel, then one can write the exact wave function ^'^ ^ as, 

« X^:nr.t)x(:nrct)- (23) 
This leads to the prior form DWBA T - matrix 

t\-;^^^^'°^\dWBA) = {x[:Hnt)xt'^{rct)\Vd + Vu - f/at|x(+Hr,,)0a(r6c)>. (24) 

It can be shown 021 that the DWBA T - matrices given by Eqs. (22) and 
(24) are equivalent to one another, i.e., 

j,j+)bo.t]p^^^^ ^yj-)[p™r]^^^^^^_ (25) 

Thus, for actual calculations, one may use the T - matrix which seems more 
convenient. y(-)[p'"«o'"] involves very complicated coordinate transformations 
as compared to the T^+^I^*"^*! form and hence it is very difficult to work with 
it in actual problems. Moreover Vt,c in Eq. (22) is of a shorter range than 
Vet + — Vat in Eq. (24), which would make the numerical evaluation of 
Eq. (22) relatively easier. The post form DWBA has been extensively used to 
perform breakup calculations [23 EH |H ESI • 

However, by introducing a different approximation for 'I'y an alternate 
prior form T - matrix can be obtained. If we assume that the final state inter- 
action Vbc between the fragments h and c is important, then one can approximate 
as 

*^^«x£^(rat)0i"!q,(r.c). (26) 
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In Eq. (26), the relative motion wave function of b and c (which could also 
be a resonant state) is described by (/^i* (rfjc), where q/ denotes the relative 
momentum between the fragments. The center of mass (cm.) motion of the 
unbound system (a* = b + c) with respect to the target in the final state is given 
by XQ^''(i"at) with momentum Qf. They are related to momenta q^ and qc of 
fragments b and c by 



and 



Q/ = qb + qc (27) 



q/ = — qb qc, (28) 



respectively. 

This approximation [Eq. (26)] leads to an alternate prior form T - matrix 

Ti-)lalt„]^^^^^^ - {X^Q^iraMi^l.incmt + Vbt - UaU;A^at)Unc)) (29) 



or 

= (xy(raO'^i"l,,(rbc)|^ct + F.*|x^+Hr<,,)^a(r.c)). (30) 

In Eq. (29), Uat depends on r^t while (/>i^^q^ (rf,c) and (jiainc) depends on r^c, 
and hence the explicit dependence on Uat in Eq- (30) has dropped out because 
of the orthogonahty of </>i~_q^ (fbc) and 4>a{Ybc)- We note that Eq. (30) contains 
only two distorted waves as against the three distorted waves in T^+)\p°^''-\ The 
distorted waves being oscillatory even at large distances, any reduction in their 
number will accelerate the convergence of the T - matrix. The alternate prior 
form has been used for performing breakup calculations by several authors |46[ 
EH. It should be noted that T^~^^''^*'^"°''\dWBA) is no longer equivalent to 

t^+)\p°^^\dwbA). 

The T - matrix [Eq. (30)] describes a situation in which the projectile a is 
inelastically excited from the ground state to its continuum. If we ignore the 
nuclear interactions in the distorted waves in both incident and final channels 
and also ignore the nuclear parts of interactions Vd and Vbt , then the alternate 
prior form T - matrix describes the Coulomb excitation of the projectile. The 
semi-classical counterpart of Eq. (30) is the Alder- Winther theory of Coulomb 
excitation 0H|- 

The alternate prior form DWBA can be regarded as the first iteration of the 
solutions of a coupled channels problem e.g., the coupled discretized continuum 
channels or CDCC equations (see, e.g. "WT). However breakup studies of both 
stable isotopes 0111201 and halo nuclei ^52 , .5311^1 have shown that the alternate 
prior form DWBA is insufficient to describe the data and that higher order 
coupling effects of the breakup channels are important in both the cases. 

In the next sub-section, we show that by suitably rewriting the post form 
DWBA T - matrix [Eq. (22)], the features of the spectator-participant mecha- 
nism can be explicitly seen therein j55| . 
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2.2 Post form DWBA T — matrix in quasi free limit 

In Eq. (22) , wc replace the relative motion wave function of the projectile by a 
plane wave, i.e., 

X^:^(rat)=e^''-"-°% (31) 

and rewrite the scattering wave functions in the final state, in terms of the half 
off-shell f-matrix elements [f^t (p, q^), £ = /2mj] as, 

xtV..) = a-- + l.n; pi)3 / -^'v-^f^^, j = 0,c. (32) 

The T - matrix then has four terms: (1) a term containing three plane waves, 
which vanishes because of energy conservation, (2) a term (T^) describing the 
scattering of particle c on t, 

^- = 7^1i°i/ / / '^VrfV„td3p^^^^l!l-4(p,q,)e-''i^---'''y6e(r6c)e*'»"-'-»*(/)a(r6c), 

(33) 

(3) a term (Tf,) [which is similar to that given by Eq. (33)] describing the 
scattering of particle b on t, (4) and a term (T(,c) describing the scattering of 
both b and c on t (double scattering). 

Let us consider, e.g., the term [Eq. (33)]. We express r^t and r(,t in terms 
of rftc and Tat (Fig. 2) as. 



b 




Figure 2: The coordinate system used in section 2.2. 
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and 



Yet = rat nc (34) 

ma 



rbt = rat H nc- (35) 

ma 



Substituting Eqs. (34) and (35) in Eq. (33) we get, 

p»(qa-q6-p)-r, 



1 r f r pi(<ia-<ib-p)-rat „^ , 

= Wrl'^oJ J jM'atd'p e-(^P-^'^-)--°t;,(p,qe) 

X ne(r6c).^a(r6c) (36) 

Integrations over coordinates rat and p can be performed easily to give, 



T, = te^^lf!£)j.(Q), (37) 

(qa-q5)^-?? 



where 



with 



F(Q) = y d''nce''^''''VUnMnc), (38) 



Q = I^q, - q,,. (39) 
Now the Schrodinger equation for the b—c system can be used to rewrite Eq. (38) 



as 

i2 



F(Q) = -{2nfn^{a''+Q^)'^am (40) 



where 



^„(Q) = (27r)-3/2 j d^bce'^^-^'-ct^aivbc) (41) 

is the momentum space wave function of the projectile in its ground state. In 
Eq. (40) a is defined as 

£60 = ^, (42) 

where € is the binding energy of the h — c system and Hf,c is the corresponding 
reduced mass. Thus is written as 

^Hbc (Qa - q&j - 
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Now, for an infinitely heavy target, the energy conservation gives 

- t£!±|!^.=i. (44) 

(qa - qtr - li "^a 

Substituting Eq. (44) in Eq. (43) we obtain, 

To = —{2TTf/XMa - qb, qc)0a(Q). (45) 

A similar expression holds for Tf,. 

If the off-shell dependence of i*t(qa — qf),qc) in Eq- (45), is neglected |55| . 
then the breakup cross section is determined by the modulus square of the wave 
function, 0a (Q)- The maximum cross section would occur when 

Q = — qa - qb = 0, (46) 

ma 

or 

Va = Vb, (47) 

i.e., when particle b moves with the beam velocity and is only a 'spectator' 
in the breakup process. Thus, the term Tc (or Th) has the simple physical 
interpretation: particle c (or b) interacts with t via Vd (or Vbt) and is knocked 
out of the projectile, while particle b (or c) is only a spectator. This shows that 
the main feature of the Serber model - the dependence of the cross section on 
the internal momentum distribution of the fragment within the projectile - is 
embedded within the post form DWBA T - matrix in a quasi free limit. 

It must, however, be cautioned that while some gross physical insights are 
obtained from the quasi free limit of the post form DWBA T - matrix, actual 
calculations must be made with full distorted waves j5()| in order to explain the 
breakup data. 

The breakup reactions of halo nuclei have been investigated by several au- 
thors using a variety of approaches (see, e.g., |57| EH EH EOl i and for an 
exhaustive bibliography on the subject). However, only a few '47'"62","63", "64II65| 
of them treat both Coulomb and nuclear breakup terms consistently on an equal 
footing. In Ref. |Ej the prior form T - matrix given by Eq. (30) has been used 
to study the breakup reactions of ^B. It should be recalled that in this T - matrix 
the fragment-target interaction is treated in first order which was subsequently 
shown to be inadequate. In Refs. jHH ES] the time evolution of the projectile 
in the coordinate space is described by solving the time dependent Schrodinger 
equation, treating the projectile target interaction as a time dependent external 
perturbation. These calculations use the semiclassical concept of the motion of 
the projectile along a trejectory. 

The post form DWBA formulation of breakup reactions which uses the T 
matrix given by Eq. (20) includes consistently both Coulomb and nuclear inter- 
actions between the projectile fragments and the target nucleus to all orders, 
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but treats the fragment-fragment interaction in first order. As can be noted 
easily, the full wave function describing the ground state structure of the pro- 
jectile, enters as an input in this theory which makes it possible to investigate 
the structure of the projectile from the study of the breakup reactions. It can 
treat the Coulomb and nuclear breakups as well as their interference terms con- 
sistently on an equal footing. Since this theory uses the post form scattering 
amplitude, the breakup contributions from the entire valence nucleon-core frag- 
ment continuum corresponding to all the multipoles and the relative angular 
momenta are included in it. This can account for the postacceleration effects 
in a unique way. In the subsequent sections we describe this theory in some 
details. 

3 Breakup amplitudes in the post form DWBA 
3.1 Pure Coulomb breakup 

Coulomb dissociation of halo nuclei has been investigated by several authors 
using a number of different theoretical approaches. A semiclassical coupled 
channels formalism has been used by authors of Ref. while in Refs. |67II68| 
the time dependent Schrodinger equation method was employed. The results 
within these approaches depend on the range of the impact parameter associated 
with the straight line trajectories used to describe the motion of the projectile 
in the field of target nuclei. However, in these studies the emphasis was on 
investigating the dynamics of the Coulomb dissociation, and not the structure 
of the projectile ground state which was assumed to have some very simple zero 
range (ZR) form. Similar assumption for the projectile structure was also made 
in other semiclassical 69. . 70; and prior form distorted wave Born approximation 
(DWBA) calculations |7T] . 

In this section, we present a theoretical model to describe the pure Coulomb 
breakup of one-neutron halo nuclei within the framework of the post form 
DWBA where finite range effects are included via a local momentum approxi- 
mation (LMA) |72l 1781 17^ . This theory of breakup reactions incorporates the 
details of the ground state structure of the projectile in the breakup ampli- 
tude [371|77|. 

We consider the reaction a + t b + c + t, where the projectile a breaks up 
into fragments b (charged) and c (uncharged) in the Coulomb field of a target 
t. The coordinate system chosen is shown in Fig. 3. 

The position vectors satisfy the following relations 

r — Ti — ari, a = — (48) 

rric + rah 

Tc = jri+Sri, (5= ■ , 7 = (l-a(5) (49) 

TOb + nit 

The starting point of our discussion is the post form DWBA T - matrix [Eq. 
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t 



Figure 3: Coordinate system used in the T-matrix (Eq. 50). b, c and t represent 
the charged core, valence neutron and target, respectively. 

(20)] 

^ = 111 rferfridrai"^*(k6,r)$^(6)x(-'*ke,re$:(^c)V6c(ri) 

X $a(ea,ri)xi+'(ka,r,), (50) 

where x's are the distorted waves for relative motions of b and c with respect to 
t and the center of mass (cm.) of the b — t system, respectively, and $'s are the 
internal state wave functions of concerned particles with internal coordinates ^. 
kb, kg are Jacobi wave vectors of b and c, respectively in the final channel of the 
reaction. Vbc{ri) is the interaction between b and c. The charged fragment b and 
the projectile a interacts with the target by a point Coulomb interaction and 
hence the corresponding distorted waves are the Coulomb distorted waves with 
appropriate boundary conditions. For pure Coulomb breakup, the interaction 
between the target and uncharged fragment c is zero and hence for this case 
we can write Xc ''*kc,rc = e^^(K-'^a) _ $(i(^a,ri) represents the bound state 
wave function of the projectile having its radial and angular parts as u^(ri) and 
F,f,(fi), respectively, which arc associated with the relative motion of b and c. 
Integrals over the internal coordinates ^ give 

= '^{imjcfJ'c\3fJ'){jbfJ'b3fJ'\3afJ'a)i%{ri)Yl{ri)- (51) 

imjfj. 

In Eq. (51), ^ is the orbital angular momentum for the relative motion between 
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b and c, ja is the spin of a. Using this equation the T-matrix can be written as 
T = {(^rnjcfic\jfJ'){jbf^bjfJ'\jaf^a)i^'iPhnOib, kc; kg), (52) 

Imjfj, 

X0r(ri)xi+nka,r.), (53) 

with (3im being the reduced T - matrix and £ = \/2i+ 1 . We have written 
/-(ri)=w,(ri)r^(fi). 

It may be noted that the reduced ampHtude (iim involves a six dimensional 
integral which makes its evaluation quite complicated. The problem gets fur- 
ther acute due to the fact that the integrand involves three scattering waves 
which have oscillatory behaviour asymptotically. Therefore, several approxi- 
mate methods have been used in the literature to avoid the evaluation of six 
dimensional integrals. In the zero range approximation (ZRA) 41 one assumes 

T4,(ri)</>^(ri) = 7^o<5(ri), (54) 

where Dq is the usual zero range constant. This approximation reduces the six 
dimensional integral in Eq. (53) to a three-dimensional one. The corresponding 
amplitude is written as 

PiJ" - Do{xt\k,,r,)e^'^-'^^\xi+H^a,r,)). (55) 

In Eq. (55), the details of the projectile structure enter in the reaction amplitude 
only as a multiplicative constant Dq. However, ZRA necessarily restricts the 
relative motion between b and c in the projectile a to s - state only. Even for 
such cases, this approximation may not be satisfied for heavier projectiles and 
at higher beam energies 

Baur and Trautmann (BT) |75| have proposed an alternative approximation 
in which the projectile cm. coordinate in the corresponding distorted wave in 
Eq. (53) is replaced by that of the core-target system, i.e. « r. With this 
approximation the amplitude splits into two factors each involving a three 
dimensional integral 

^>fJ = (e''^-'-MT4c|0r(ri))(xr^(k,,r)e^^''-nxi+Hka,r)). (56) 

The first term depends on the structure of the projectile through its ground 
state wave function 0^™(ri). The second term involves the dynamics of the 
reaction, which can be expressed analytically in terms of the bremsstrahlung 
integral |S7]. This amplitude (which will be referred to as the BT amplitude), 
used originally to study the deuteron breakup at sub-Coulomb energies |75| . 
was applied to calculations of the Coulomb breakup of halo nuclei in Ref. |76|. 
This approximation, which allows the application of the theory to non-s - wave 



where 



3£m(kf„kc;ka) ^ J J 
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projectiles, may seem to be justified if the cm of the b — c system is shifted 
towards b (which is indeed the case if nib ^ rric). However, the neglected piece 
of Ti (ari ) occurs in association with the wave vector , whose magnitude could 
be quite large for the reactions at higher beam energies. Therefore, contributions 
coming to the amplitude from the neglected part of may still be substantial. 

An approximate way of taking into account the finite range effects in the 
post form DWBA theory is provided by the local momentum approximation 
| 74l I72| . The attractive feature of this approximation is that it leads to the 
factorization of the amplitude Pim similar to that obtained in the BT case. We 
use this approximation to write the Coulomb distorted wave of particle b in the 
final channel as 

xi~\^b,r) = e-^"K-"-^x^^(k,,r,). (57) 

Eq. (57) represents an exact Taylor series expansion about if K(= — iVrJ 
is treated exactly. However, this is not done in the LMA scheme. Instead, the 
magnitude of the local momentum is taken to be 

K{R) = ^^[E-ViR)], (58) 

where m is the reduced mass of the b — t system, E is the energy of particle b 
relative to the target in the cm. system and V{R) is the Coulomb potential 
between b and the target at a distance R. Thus, the local momentum K is 
evaluated at some distance R, and its magnitude is held fixed for all the values 
of r. As shown in appendix of Ref. [ZJ], the magnitude of K remains constant 
for r > 10 fm. Due to the peripheral nature of breakup reactions, this region 
contributes maximum to the cross section. Therefore, we have taken a constant 
magnitude for K evaluated at i? = 10 fm for all the values of r. As is discussed 
in [231^ the results of our calculations are almost independent of the choice of 
the direction of the local momentum. Therefore, we have taken the directions 
of K and kb to be the same in all the calculations. Detailed discussion on the 
validity of the local momentum approximation is presented in Refs. [771 164| . 
It may be noted that in the calculations presented in Ref. [TSj, the LMA was 
applied to the Coulomb distorted wave of the projectile channel which could 
imply a deviation from the distorted wave approximation. 

Substituting Eq. (57) into Eq. (53) we get the following factorizcd form of 
the reduced amplitude 

= (e'('''^-"^)-'-M^fcc|</'i"(ri)) 

X (x^^(k,,r,)e'*''-"--|xi+Hk„r,)). (59) 

Eq. (59) (which will be referred to as the FRDWBA amplitude in the following) 
looks like Eq. (56) of the BT theory but with the very important difference that 
the form factor is now evaluated at the momentum transfer (7kc — aK), and not 
at the valence particle momentum kc. The two momenta could be quite different 
for cases of interest in this work. The second term, involving the dynamics of 



15 



the reaction, is the same in both the cases. Therefore, the breakup amphtude 
obtained in BT approximation differs from that of FRDWBA by a factor 

p ^ ^ (e-'--'-M^6cie'(ri)) 




Figure 4: Modulus square of the ratio Fr defined by Eq. (60) as a function of 
the beam energy for the breakup of ^^Be on a gold target at the beam energy 
of 44 MeV/nucleon corresponding to forward emission angles of the breakup 
fragments as discussed in the text. 

In Fig. 4, we have shown the beam energy dependence of for the 

breakup of ^^Be on a gold target at the beam energy of 44 MeV/nucleon for a 
set of forward angles of the outgoing fragments {9i, — 1°, 6c = 1° and (f>c = 1°). 
In this calculation, the ground state wave function of ^^Be has been obtained by 
considering a configurations in which as- wave valence neutron is coupled to 0^ 
^°Be core [^°Be(0+) (g) lsi/2i^] with the one-neutron separation energy (Sn-wg^) 
of 504 keV and a spectroscopic factor (SF) of 0.74 [H EOl EI- The single 
particle wave function has been obtained by assuming the valence neutron-^'^Be 
interaction to be of Woods-Saxon type whose depth is adjusted to reproduce 
the corresponding binding energy with fixed values of the radius and diffuseness 
parameters (taken to be 1.15 fm and 0.5 fm, respectively) (see [77] for more 
details). In this figure, we see that |i^rP is close to unity only at the sub- 
Coulomb beam energies (of course at higher beam energies it crosses twice the 
line representing the value 1). Therefore, the BT and FRDWBA calculations 
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are expected to produce similar results at very low incident energies. Depending 
upon the beam energy, the BT results can be larger or smaller than those of the 
FRDWBA. 

Recently, a theory of the Coulomb breakup has been developed within an 
adiabatic (AD) model |S21IHSli where one assumes that the excitation of the 
projectile is such that the relative energy (ebc) of the b — c system is much 
smaller than the total incident energy, which allows ebc to be replaced by the 
constant separation energy of the fragments in the projectile ground state. It 
was shown in [53] that under these conditions the wave function \l/i^^(^a, ri, rO 
has an exact solution as given below 

*i-*'^'^''(Ca,ri,r,) = a>,(ea,ri)e^"''-'-^xi+Hk„r). (61) 
Substituting 'J'i^''''^^ for 'J'i^'' in Eq. (50) leads to the reduced amplitude: 
if3if = (e'('^-"'^»)-Mnc|/r(ri))(x^Hk,,r)e^^''-nxi+Hk<„r)). (62) 

This amplitude differs from those of BT as well as FRDWBA only in the form 
factor part (first term), which is evaluated here at the momentum transfer 
(kc — aka). Eq. (62) can also be obtained in the DWBA model by making 
a local momentum approximation to the Coulomb distorted wave in the initial 
channel of the reaction, and by evaluating the local momentum at i? — oo with 
its direction being the same as that of the projectile 78 . The adiabatic model 
does not make the weak coupling approximation of the DWBA. However, it 
necessarily requires one of the fragments (in this case c) to be chargeless. In 
contrast, the FRDWBA can be applied to cases where both the fragments b 
and c are charged (see, e.g., Ref. [H]). While the effect of nuclear breakup in 
the adiabatic model description of the elastic scattering of the loosely bound 
projectile has been calculated in Refs. ,83-'_84.._85,7 the nuclear part of the ampli- 
tude for breakup reactions is yet to be calculated within this model. However, 
the calculations of the nuclear breakup cross section has been done within the 
FRDWBA theory 64 . 

The triple differential cross section of the reaction is given by 

'^'^ ^p{Eb,nb,^lc)y\Ptm\^ (63) 



dEbdftbdilc tiva ' ' 

where p{Eb, fib, ^Ic) is the appropriate [SS] three-body phase space factor, given 

by 

n n \ h-^mbnicmtPbPc 

mt+nic- nic — ^-p — — 

with ka,kt, and kc being evaluated in the appropriate frame of reference. Va 
is the relative velocity of the projectile in the initial channel. In Eq. (64), the 
linear momenta, p, are related to wave numbers k hy p = fik. 
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Substituting the following expressions for the Coulomb distorted waves 
xi"^*(kb, r,) = "'/2r(l + i776)e-''''-'-'iFi(-«yb, l,i{hn + kfc.r0),(65) 

X^+\K,r^) = e-^^-/^r{l + iT^a)e'''-''\Fi{-irja,l,t{kaU-ka.n)) (66) 
in Eq. (59), one gets for the triple differential cross section 



dEhdrtbdric fiVa 



In Eqs. (65 ~ 67), 77's are the Coulomb parameters for the concerned particles. 
In Eq. (67), / is the bremsstrahlung integral [ST] which can be evaluated in a 
closed form: 



-I 

'dB 



(-) 

\ dx / x=Q 



dx ) x=o 

iFi{-i7ja,-iVb;l;D{0)) 



(68) 



where 

B{x) 



A-K 



(fc2 - 2k.ka - 2xfca)"'°(fc^ - 2k.k6 - 2xfc6)*'"' , (69) 



D{x) 



2P{kakb + kg.kfc) - 4(k.ka + Xka){k.kb + Xkb) 

(fc2 - 2k.ka - 2xka){k^ - 2k.kfc - 2xkb) 



(70) 



with k = ka — kf, — (5kc. Zi contains the projectile structure information and is 
given by 



Zi^ drirlji{kiri)Vbciri)uiiri), 



(71) 



with ki = |7kc — aK|. 

Let us now discuss some numerical applications of the theory of the pure 
Coulomb breakup reactions as presented above. We investigate the breakup of 
neutron rich nuclei ^^Be and i5.i'i',i9(^ j^g^]^ energies below 100 MeV/nucleon. 
Apart from the distance at which the local momentum is calculated (which is 
taken to be 10 fm ) and its direction (described earher), the only other input to 
our calculations is the radial part of the projectile ground state wave function. 
As discussed above, we have assumed a Woods-Saxon potential to describe the 
valence neutron-core relative motion whose depth is searched, for a given con- 
figuration, to reproduce the corresponding binding energy. In the calculations 
presented here we have mostly considered a s-wave configuration (as described 
above) for the ^^Be ground state. However, in some cases we have also consid- 
ered a c?-wave configuration for this nucleus in which a rf-wave valence neutron 
is coupled to 2+ ^°Be core [^'^Be(2+) (g) 6^5/2^ with the one- neutron separation 
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energy {Sn-^OBe) of 3.872 MeV [77]. The configurations for the C isotopes used 
in our calculations are described at the appropriate places below. 

In Fig. 5, we present a comparison of our calculation with the data (taken 
from [7D]) for the neutron energy distribution of the double differential cross 
section {(Pa/dEnd^ln) at the neutron angles of 1° and 3.4°, in the breakup of 
^^Be on a gold target at the beam energy of 44 MeV/nucleon. Calculations 
performed within both FRDWBA and AD model of pure Coulomb breakup 
are shown in this figure. The s-wave configuration has been used for the ^^Be 
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Figure 5: Neutron energy distributions for the breakup of ^^Be on Au at beam 
energies of 37 MeV/nucleon (dotted lines) and 44 MeV/nucleon (solid lines), 
calculated using configuration (a) with single particle wave functions within the 
FRDWBA and the AD models. The top half of the figure is for 6i„ = 1°, while 
the bottom half is for 0„ = 3.4°. The experimental data are taken from |7U| . 

ground state in both the cases. The beam energy in this experiment |7()| varies 
between 36.9 - 44.1 MeV/nucleon. To take into account this spread, we have 
performed calculations at both its upper (44 MeV/nucleon) (solid line) and 
lower ends (37 MeV/nucleon) (dotted line). Even though these data have large 
statistical errors, the calculations performed at 44 MeV/nucleon are in better 
agreement with the experimental values. It should also be noted that the AD 
model calculations over-predict the experimental cross sections in the peak re- 
gion. 

The measured neutron angular distribution in the exclusive ^-'^Be + A — > 
""^^Be + n +A reaction on heavy targets below the grazing angle is very narrow 
and is shown to be |78l IHHj dominated by the Coulomb breakup process. This 
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reflects the narrow width of the transverse momentum distribution of the valence 
neutron in the ground state of ^^Be, which is consistent with the presence of 
a neutron halo structure in ^^Be. In Fig. 6, we compare the calculated and 
measured exclusive neutron angular distribution da/dfln as a function of the 
neutron angle 0„ for the above reaction on Au, Ti and Be targets at the beam 
energy of 41 MeV/nucleon. The ^^Be ground state wave function is the same 
as described above. We note that for Au and Ti targets pure coulomb breakup 
calculations are able to describe the data at forward angles (for neutron angles 
below 25° and 15°, respectively in the two case). On the other hand, for the 
^Be target pure breakup calculations are much below the data everywhere. This 
gives a clear indication of the importance of nuclear breakup effects at the 
backward angles for the medium mass and heavy targets and everywhere for 
the light target. 
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Figure 6: The calculated neutron angular distributions for the breakup of ^^Be 
on a Au, Ti and Be targets at 41 MeV/nuclcon beam energy. The data are 
taken from |70| . 



Pure Coulomb contribution to the relative energy spectrum in the breakup 
of ^^Be on a Pb target at 72 MeV/nucleon is shown in Fig. 7. The top half 
shows the results obtained with FRDWBA (solid line), AD model (dotted line) 
and the BT approximation (dashed line) using the s-wave configuration for 
the iiRe ground state. In the bottom half, we show additionally the results 
obtained with the ^^^Be wave function calculated within a dynamical core po- 
larization (DCP) model [HHI (dotted line) and that obtained with the c?-wave 
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configuration as described above. We see that, while both the FRDWBA and 
the AD model calculations reproduce the peak value of the data |90j well, the 
FRDWBA calculations done with the DCP wave function overestimate it. On 
the other hand, none of the calculations is able to explain the data at higher 
relative energies. This can be attributed to the fact that nuclear breakup effects, 
which can contribute substantially at higher relative energies (for Erei > 
0.6 MeV), are not included in these calculations. Of course, authors of Ref. 
PU) claim that their data have been corrected for these contributions. However, 
the procedure adopted by them for this purpose is inadequate. They obtained 
the nuclear breakup contribution on the Pb target, by scaling the cross sections 
measured on a carbon target. This scaling procedure is unlikely to be accurate 
for reactions induced by halo nuclei due to the presence of a long tail in their 
ground state. 
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Figure 7: Relative energy spectra for the Coulomb breakup of ^^Be on a Pb 
target at 72 MeV/nucleon beam energy. The top half of the figure shows the 
spectra obtained with a single particle wave function, using the FRDWBA (solid 
line) , the AD model (dotted line) and the BT approximation (dashed line) . The 
bottom half shows the results of FRDWBA calculations using single particle 
(solid line) and DCP (dotted line) wave functions. The dot-dashed line shows 
the d - state FRDWBA calculation using a single particle wave function, after 
being multiplied by a factor of 5000. The data are taken from |9()| . 

In a full quantum mechanical theory, both Coulomb and nuclear breakup 
contributions should be calculated on the same footing and corresponding am- 
plitudes should be added coherently. This is discussed in the next section. 
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Calculations done using the BT theory (dashed line in the upper part of Fig. 7) 
underestimates the data considerably. This difference between the FRDWBA 
and the BT results can again be traced to the behaviour of l-FrP in Fig. 4, which 
is smaller than unity at the beam energy of 72 MeV/nucleon of this reaction. 
The FRDWBA result using the d - state configuration with single particle wave 
function (dot-dashed line in the lower half of Fig. 7, shown after multiplying 
the actual numbers by 5000) also grossly underestimates the data. 

The neutron halo structure is reflected in the narrow width of the parallel 
momentum distribution (PMD) of the charged breakup fragments emitted in 
breakup reactions induced by the halo nuclei. This is because the PMD is least 
affected by the reaction mechanism |^ |SI 1^^! IHSI and therefore, a narrow 
PMD can be related to a long tail in the matter distribution in the coordinate 
space via Heisenberg's uncertainty principle. 




250 



(Mev/c) 

Figure 8: Parallel momentum distributions of ^°Be in the breakup of ^^Be on 
Ta (top half) and U (bottom half) at the beam energy of 63 MeV/nucleon, in 
the rest frame of the projectile, with s-wave (solid line) and d-wave (dot-dashed 
line) configurations for the ^'^Be ground state. Both calculations are done in 
FRDWBA and are normalized to the data peaks. The data are taken from I93| . 

In Fig. 8, we present the PMD of the ^°Be fragment emitted in the breakup 
of ^-'^Be on U and Ta targets at 63 MeV/nucleon beam energy. Calculations 
performed within the FRDWBA and the AD models using both the s-wave 
and d-wave configurations are presented in this figure. The calculated cross 
sections are normalized to match the peak of the data points (which are given 
in arbitrary units) the normalization constant being the same for both 

cases. With the s-wave configuration the full width at half maximum (FWHM) 
for the U and Ta targets are 44 MeV/c and 43 MeV/c, respectively in both the 



22 



FRDWBA and the AD cases. These agree well with the averaged experimental 
value of 43.6±1.1 MeV/c |M| and also with those calculated in [7S|. The very 
narrow widths of the parallel momentum distributions signal the presence of a 
neutron halo structure in ^^Be. It may be noted that the PMD calculated with 
a pure d - wave configuration is too wide in width and grossly overestimates 
the experimental FWHM. 

In Fig. 9, we present the PMD (calculated within the FRDWBA formalism) 
of the ^^C fragment in the breakup of ^^C on a Ta target at the beam energy 
of 88 MeV/nucleon. We have normalized the peaks of the calculated PMDs to 
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Figure 9: FRDWBA results for the parallel momentum distribution of in 
the breakup of ^^C on Ta target at the beam energy of 88 MeV/nucleon. The 
top half shows the results obtained with the configuration [^^C (O"*") (g) IS1/2H 
and single particle wave function for the ground state of "'^^C with one-neutron 
separation energies of 530 keV (solid line), 160 keV (dashed hue). The dotted 
line shows the result obtained with a DCP wave function with a one-neutron 
separation energy of 160 keV. The bottom half shows the result obtained with 
the configurations [^^€(0+) «) IS1/2H (solid) and [i^C(0+) «) 0^5/2^ (dashed), 
with the same value of the one- neutron separation energy (530 keV). The data 
have been taken from i94l| . 

that of the data (given in arbitrary units) [H^ (this also involves coinciding the 
position of maxima of the calculated and experimental PMDs). As can be seen 
from the upper part of this figure, the experimental data clearly favour S^^isq 
= 0.53 MeV with the s - wave n-^^C relative motion in the ground state of ^^C 
(solid line). The results obtained with the s - wave configuration within the 
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simple potential (dashed line) and DCP (dotted line) models (with the same 
value of Sn-wc = 0-16 MeV) are similar to each other. 

In the lower part of Fig. 9, we have shown the results obtained with the d - 
wave relative motion (dashed line) for this system (with Sn-^>ic = 0-53 MeV) 
and have compared it with that obtained with as- wave relative motion (solid 
line) with the same value of the binding energy. As can be seen, the FWHM of 
the experimental PMD is grossly overestimated by the d - wave configuration. 
The calculated FWHM with the s - state configuration (with Sn-isc = 530 
keV) is 40 MeV/c, which is in excellent agreement with the experimental value 
of 41±3 MeV/c Thus these data favour a configuration [^^€(0+) » lsi/2i'], 
with a one-neutron separation energy of 0.530 MeV for the ground state of ^^C. 
These results are in agreement with those of Ref . . The narrow width of the 
PMD provides support to the presence of a one-neutron halo structure in ^^C. 

We next consider the breakup of ^^C which has a relatively larger value for 
the one-neutron separation energy (1.2181 MeV) and a ground state spin-parity 
of 1/2+ [SJ- We consider two configurations: a lsi/2 neutron coupled to a ^''C 
(0+) core and a 0^5/2 neutron coupled to a ^"^C (0+) core. One could have also 
considered a ^'^C (2+) core and 0^5/2 neutron coupling to get a 1/2+ ground 
state for ^^C, but it would raise the one-neutron separation energy to about 7.01 
MeV, which is highly unfavourable for the formation of a halo. We, therefore, 
do not consider this configuration in our study. 
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Figure 10: Parallel momentum distributions of ^"'C in the breakup of ^^C on Ta 
at 84 MeV/nucleon. The solid line and dotted lines show the results obtained 
with the configurations [^"'C (0+) (g) lsi/2i^] and [^''C (0+) (8) 0^5/21^], respectively 
for the ground state of the projectile. The data have been taken from |94j . 
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Table 1: FWHMs from the parallel momentum distribution of for different 
ground state configurations of ^^C and one-neutron separation energies (e) in 
the breakup of ^^C on Ta at 84 MeV/nucleon beam energy. 



Projectile 


e 


FWHM 


config 


juration 


(MeV) 


(MeV/c) 


^''C (0+) 


(g) \Sxl2V 


0.729 


51 


(0+) 




0.729 


114 


(2+) 




2.5 


82 


16c (2+) 




2.5 


185 



In Fig. 10, we present the results of our calculations for the PMD of the 
i^'C fragment in the breakup of i^C on a Ta target at the beam energy of 
84 MeV/nucleon. The experimental data are taken from [HI]. The s - state 
configuration for the ground state of i^C gives a FWHM of 62 MeV/c, while 
with the d - state configuration it comes out to be 140 MeV/c. Therefore, 
the experimental value for the FWHM (67±1 MeV/c) [M] favours the former 
configuration. Thus, our results provide support to the existence of a halo 
structure in i^C. This nucleus provides an example of the one-neutron halo 
system with the largest one-neutron separation energy, known so far. 

I'^C has a lower one-neutron separation energy (729 keV) as compared to 
that of i^C. It would be interesting to see if it also has a halo structure, which 
seems probable if one considers only the binding energies. The quoted ground 
state spin-parities for this nucleus are l/2+,3/2+ and 5/2+ [HHI- RMF cal- 
culations [nni predict it to have a value of 3/2+. We consider four possible 
ground state configurations for this nucleus and calculate the parallel momen- 
tum distributions of the i^C fragment in the breakup of i^C on a Ta target at 
84 MeV/nucleon beam energy within our FRDWBA formalism. The FWHMs 
of the PMD obtained with different configurations are listed in Table 1. 

It is evident from this table that the s - state configurations predict a narrow 
width for the PMD, providing support to the existence of a halo structure in 
this nucleus. The experimental data |H2 for the breakup of i^C on a light target 
(Be) at 84 MeV/nucleon give a FWHM of 145±5 MeV/c. Since the width of the 
PMD is mostly unaffected by the reaction mechanism , it is quite likely that 
the experimental FWHM will be the same also for the breakup of this nucleus 
on a heavier target. Therefore, the results shown in Table 1 seem to provide 
support to a d - wave configuration for the ground state of i^C |100) . Hence, 
the existence of a one-neutron halo structure is quite improbable in i^C (see 
alsoRefs. |TTrn[TT7^ V 

In summary, from our pure Coulomb breakup studies we can say that the nu- 
clei iiBe, i^C and i^C have a one-neutron halo structure in their ground states. 
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However, for ^'^C such a structure appears to be less likely. The calculated cross 
sections are selective about the ground state wave function of the projectile. At 
the same time there is a clear indication of the fact that only pure Coulomb or 
pure nuclear breakup calculations may not be sufficient to describe the details 
of the halo breakup data. Consideration of both modes of breakup as well as of 
their interference terms is necessary to describe the data properly. 

3.2 Full breakup amplitude 

The full breakup amplitude that includes consistently both Coulomb and nuclear 
interactions between the projectile fragments and the target nucleus to all orders 
has been developed in Refs. [^111^. We perform a Taylor series expansion of 
the distorted waves of particles b and c about and write 

xt\k,,r) - e-'^"K-'-^xi~^(kb,r,), (72) 
x(-)(kc,r.) = e^^K-'-^xi-Hkc,<5r,). (73) 

Employing the LMA |74[ 17^ , the magnitudes of momenta Kj are taken as 



K,{R) = ^{2m,/h^)[E,-V,(R)], (74) 

where nij {j = b, c) is the reduced mass of the j — t system, Ej is the energy of 
particle j relative to the target in the center of mass (cm.) system, and Vj{R) 
is the potential between j and t at a distance R. Substituting Eqs. (72) and 
(73) in Eq. (50), and introducing the partial wave expansion of the distorted 
waves and carrying out the angular momentum algebra, one gets 



X yL'ikb,h){LbOL,0\LaO)nL,,L^,LAka,h,ka), (75) 



where 



ytlihX) = Y.^-f'^^bML,~M\LaQ)YMikb)Y^f*{k), (76) 

M 

x: 

Zi{Q) = I rldnje{Qn)ue{n)Vbc{ri), (77) 



f°° rlr 

/ —fLAka,r,)fLAkb,r,)fLAkc,Sr,). (78) 
Jo 



In Eq. (75), Q is the magnitude of the vector Q — — aKf,. Functions / 
appearing in the radial integral 'R-Li,,La,La{ka, fcb, kg) are the radial parts of the 
distorted wave functions x's- These are obtained by solving the Schrodinger 
equation with appropriate optical potential which include both Coulomb and 
nuclear terms. The slowly converging radial integral 'R-Lb,L^,L^ are effectively 
handled by using the complex plane method |1(J3[ 111)4] . 
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This theory can be used to calculate breakup of both neutron and proton 
halo nuclei. Generally, the maximum value of the partial waves La, Lb, Lc must 
be very large in order to ensure the convergence of the partial wave summations 
in Eq. (75). However, for the case of one-neutron halo nuclei, one can make use 
of the following method to include summations over infinite number of partial 
waves. We write Pim as 

L^^^ oo 

f3em = ML,) + (79) 

where /3 is defined in the same way as Eq. (75) except for the summation sign 
and Li corresponds to Lq, Lb, and Lc- If the value of is chosen to be 

appropriately large, the contribution of the nuclear field to the second term of 
Eq. (79) can be neglected and we can write 

oo 

Plm{Li) 

where the first term on the right hand side, is the pure Coulomb breakup ampli- 
tude which for the case where one of the outgoing fragments is uncharged, can 
be expressed analytically in terms of the bremsstrahlung integral (see, e.g., Ref. 

Therefore, only two terms, with reasonable upper limits, are required to 
be evaluated by the partial wave expansion in Eq. (79). 

In the numerical applications for the breakup of ^^Be, the structure func- 
tion Zg has been calculated with the s-wave configuration for the ^^Be ground 
state. The neutron-target optical potentials were extracted from the global set 
of Bechhetti-Greenlees (see, e.g, |1U5) V while those for the ^°Be-target system 
were taken from f'105l [TnS| ). Following we have used the sum of these two 
potentials for the ^^Be-target channel (see, [BJ for more details). We found 
that values of Lf*"^ of 500 for Au, Ta, U, Pb and Ti targets and 150 for Be and 
C targets provided very good convergence of the corresponding partial wave ex- 
pansion series. The local momentum wave vectors are evaluated at a distance, 
i? = 10 fm, in all the cases with their directions being the same as that of 
asymptotic momenta. 

In Fig. 11, we show our results for the neutron angular distributions (da/dfln) 
for the breakup of ^^Be on Au, Ti and Be targets at 41 MeV/nucleon. The 
neutron energy has been integrated from 26 MeV to 80 MeV, and the core scat- 
tering angle in the lab system {9b) has been integrated from 0° to 30° for the 
Au target case and from 0° to 20° for Ti and Be target cases. The dotted and 
dashed lines represent the pure Coulomb and nuclear contributions, respectively 
while their coherent sums are shown by the solid lines. The plus signs and the 
inverted solid triangles represent the magnitudes of the positive and negative 
interference terms, respectively. Our calculations are in good agreement with 
the experimental data (shown by solid circles) for all the three targets. 

For the Be target, da /(Kin is governed solely by the nuclear breakup effects 
at all the angles. The pure Coulomb breakup contributions are down by at least 
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an order of magnitude at the forward angles and by 2-3 orders of magnitude at 
the backward angles. The CNI terms are also small in this case. 

On the other hand, for Ti and Au targets the Coulomb terms are dominant 
at the forward angles while the nuclear breakup effects are important at larger 
angles. This is to be expected for high-Z targets, as the strong Coulomb field 
causes the fragile halo system to breakup at large distances (and hence large im- 
pact parameters), leading to the predominance of Coulomb breakup at forward 
angles. The nuclear breakup assumes importance when the breakup occurs near 
the target nucleus, consequently leading to large scattering angles. 
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Figure 11: Neutron angular distribution for the breakup reaction ^^Be + A — > 
""^^Be -f n + A at the beam energy of 41 MeV/nucleon. The dotted and dashed 
lines represent the pure Coulomb and nuclear contributions, respectively while 
their coherent sums are shown by the solid lines. The plus signs and the inverted 
solid triangles represent the magnitudes of the positive and negative interference 
terms, respectively. The data are taken from |7(Jj . 
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Magnitudes of the CNI terms vary with angle; for many forward angles they 
almost coincide with those of the nuclear breakup while at the backward angles 
they are closer to the pure Coulomb breakup contributions. Signs of these terms 
also change with the neutron angle; a feature common to all the three targets. 
It is clear that the interference terms are not negligible for Ti and Au targets 
at the forward angles. For On < 10°, the magnitudes of the CNI contributions 
are similar to those of the pure nuclear terms, leading to a better description of 
the data in this region. 
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Figure 12: The differential cross section as a function of the relative energy 
of the fragments (neutron and ^"Be) in the breakup reaction of ^^Bc on ^'^^Pb, 
^*Ti and ^^C targets at 72 MeV/nucleon. The dotted and dashed lines represent 
the pure Coulomb and nuclear breakup contributions, respectively while their 
coherent and incoherent sums are shown by the solid and dot-dashed lines, 
respectively. The plus signs and the inverted triangles represent the magnitudes 
of the positive and negative interference terms, respectively. The data are taken 
from pi) . 

The relative energy spectrum of the fragments (neutron and ^"Be) emitted 
in the breakup of ^^Be on ^osp]-, ^^^^ panel), '^^Ti (middle panel) and ^^C (bot- 
tom panel) targets at the beam energy of 72 MeV/nucleon is shown in Fig. 12. 
In these calculations the integration over the projectile cm. angle (^,iioBo~Pb) 
has been done in the range of 0°~40°, mainly to include the effects of nuclear 
breakup coming from small impact parameters. The relative angle between 
the fragments (^„_iobo) has been integrated from 0° to 180°. The dotted and 
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dashed lines represent the pure Coulomb and nuclear breakup contributions, 
respectively while their coherent and incoherent sums are shown by the solid 
and dot-dashed lines, respectively. The plus signs and the inverted triangles 
represent the magnitudes of the positive and negative interference terms, re- 
spectively. 

In case of breakup on a heavy target (^°*Pb) [Fig. (12) (top panel)] the pure 
Coulomb contributions dominate the cross sections around the peak value, while 
at larger relative energies nuclear breakup gains importance. This is attributed 
to the different energy dependence of the two contributions [HH. The nuclear 
breakup occurs when the projectile and the target nuclei are close to each other. 
Its magnitude, which is determined mostly by the geometrical conditions, has 
a weak dependence on the relative energy of the outgoing fragments beyond a 
certain minimum value. In contrast, the Coulomb breakup contribution has a 
long range and it shows a strong energy dependence. The number of virtual 
photons increases for small excitation energies and hence the cross sections rise 
sharply at low excitation energies. After a certain value of this energy the cross 
sections decrease due to setting in of the adiabatic cut-off. The coherent sum 
of the Coulomb and nuclear contributions provides a good overall description of 
the experimental data. The nuclear and the CNI terms are necessary to explain 
the data at larger relative energies. 

In the middle panel of Fig. (12), we show the relative energy of the fragments 
in the breakup of Be on a medium mass target ("''^Ti). At low relative energies 
the pure Coulomb contributions are slightly higher than the pure nuclear ones, 
while at higher relative energies it is the nuclear part which dominates. Apart 
from the very low relative energy region the CNI terms play an important role, 
which is clearly borne out by the difference in the coherent (solid) and incoherent 
(dot-dashed) sums of the pure Coulomb and pure nuclear contributions. 

The relative energy spectra for the breakup on a light target (^^C) is shown 
in the bottom panel of Fig. (12). In this case we have used the same optical 
potential for the ^°Be-^^C system as in the ^°Be-^Be case, which we had used 
earlier in calculating the neutron angular distribution in Ref. |63| . The total 
cross section in this case is normalized to the experimental cross section (found 
by integrating the area under the data points) and the same normalization 
constant is used for all the cross sections in this case. The breakup is clearly 
seen to be nuclear dominated at all relative energies, and the pure Coulomb and 
CNI terms have very little contributions. 

The parallel momentum distributions (PMDs) of the ^°Be fragment in the 
breakup of ^^Be on U and Ta targets, at 63 MeV/nucleon beam energy are 
presented in the rest frame of the projectile, in Fig. 13. The dotted and dashed 
lines show the contributions of the pure Coulomb and nuclear breakups, respec- 
tively, while their coherent sums are represented by solid lines. The coherent 
sum is normalized to the peak of the data, which are given in arbitrary units, 
and the same normalization factor has been used for the pure Coulomb and 
pure nuclear contributions. 

The above results make it clear that the Coulomb-nuclear interference terms 
are both energy and angle dependent. They are almost of the same magnitude 
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Figure 13: The parallel momentum distribution of the core in the breakup of 
^^Be on U and Ta targets, at 63 MeV/nucleon beam energy, in the rest frame 
of the projectile. The dotted and dashed lines represent the pure Coulomb and 
nuclear breakup contributions, respectively while their coherent sums are shown 
by solid lines. The data are taken from |l()7j . 



as the nuclear breakup contributions in the neutron angular distributions on 
heavy and medium mass targets. This leads to a difference in the coherent 
and incoherent sums of the Coulomb and nuclear terms, more so at forward 
angles. In the parallel momentum distribution of the ^"Be fragment in the 
breakup reaction of ^^Be the region around the peak of the distribution, which 
gets substantial contributions from forward scattered fragments, is Coulomb 
dominated, while in the wings of the distribution, where contributions come 
from fragments scattered at large angles, the nuclear breakup contributions 
dominate. 

The relative energy spectra of the fragments (neutron and ^'^Be) are largely 
nuclear dominated for light targets. However, to explain satisfactorily the data 
on heavier targets one requires both the nuclear and the CNl terms particularly 
at higher relative energies. In case of breakup on a medium mass target, the total 
pure Coulomb and pure nuclear contributions were nearly equal in magnitude. 
Thus in many sophisticated experiments planned in the future one has to look 
into the role played by the CNI terms in analyzing the experimental data. 
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4 Postacceleration effects in the Coulomb breakup 
of neutron halo nuclei 



As shown section 3.1, an important advantage of the post form DWBA the- 
ory of breakup reactions is that it can be solved analyticaUy for the case of 
the breakup of the neutron halo nuclei with the entrance and outgoing chan- 
nels involving only the Coulomb distortions |1U8[ 175] . It constitutes an ideal 
"theoretical laboratory" to investigate the physics of the breakup reactions, its 
certain limiting cases, and its relation to other models like the semiclassical 
approximation. Particularly, the effect of postacceleration can be studied in a 
unique way within this approach. 

Postacceleration refers to the situation where the core c has a larger final 
state energy than what one gets from sharing the kinetic energy among the 
fragments according to their mass ratio. This effect arises in a purely clas- 
sical picture 109 of the breakup process. The nucleus a = {c + n) moves 
up the Coulomb potential, loosing the appropriate amount of kinetic energy. 
At an assumed "breakup point" , this kinetic energy (minus the binding en- 
ergy) is supposed to be shared among the fragments according to their mass 
ratio (assuming that the velocities of c and n are equal). Running down the 
Coulomb barrier, the charged particle c alone (and not the neutron) gains back 
the Coulomb energy, resulting in its postacceleration. Of course this picture is 
based on the purely classical interpretation of this process, and will be modified 
in a quantal treatment, where such a "breakup point" does not exist. Postac- 
celeration is clearly observed in the low energy dcuteron breakup, both in the 
theoretical calculations as well as in the corresponding experiments (see, e.g., 
|37l IllUj l. However, in the description of the Coulomb dissociation of halo nu- 
clei at high beam energies within this theory [771 176L llllj , the postacceleration 
effects become negligibly small. We shall investigate this point further for the 
^^Be and ^^C Coulomb dissociation experiments '90' 112'. On the other hand, in 
the semiclassical Coulomb excitation theory the higher order effects have been 
found |113| to be small, for both zero range as well as finite range wave functions 
of the c + n system. 

It was recently noticed 61 that in the limit of Coulomb parameter <C 1 
(i.e. in the Born approximation), the pure Coulomb post form DWBA [Eq. 
(53)] leads to results which are same as those obtained in a semiclassical model 
[Tl4j . This agreement is also valid for arbitrary values of rja a-nd rjc , provided the 
beam energies are high as compared to the relative energy (Ecn) of fragments 
c and n in the ground state of the projectile. The first order approximation to 
the amplitude given by Eq. (53), can be written as [Il4. .115. .Gl, 



^ofirst order 




1 



1 



X 



2 



32 



11l 



Be + Pb 
72 MeV/nucleon 




12 







30 MeV/nucleon 



ec=en=3 




700 710 720 730 740 280 290 300 310 320 



2 
1.6 



1 MeV/nucleon 




0.4 
0.3 
0.2 
0.1 




5 MeV/nucleon 



e,=e,=25^ 




90 95 100 105 110 44 46 48 50 52 54 
E, (MeV) E, (MeV) 



Figure 14: Triple differential cross section as a function of the energy of Be 
core for the reaction ^^Be + Pb ^ n + ^°Be + Pb at the beam energies of 
72 MeV/nucleon, 30 MeV/nucleon, 10 MeV/nucleon and 5 MeV/nucleon. The 
results of the finite range DWBA and first-order theory are shown by solid and 
dotted lines respectively. 



where the relative phase (j) — cr{r]c) — cr{r]a) — cr(f) — ^/[21og |D(0)|], with air]) 
being the usual Coulomb phase shifts, and ^ = t]c — ija and D{0) as defined 
in Eq. (70). In Eq. (81), we have defined /cquI ~ "^Vaka/k^- This term is 
very similar to the Born approximation (BA) result given in |116) : in the limit 
^ ^ it actually coincides with the BA expression. This equation can be used 
to investigate the role of higher order effects (which includes postacceleration) . 
It may be noted that the derivation of Eq. (81) makes use only of the condition 
—D{0) 3> 1 which is met for beam energies large as compared to the binding 
energy. 

We now investigate postacceleration effects in the breakup of the one- neutron 
halo nuclei ^^Be and ^^C. We take a heavy target of atomic number Z=82. In 
the following all the higher order results correspond to calculations performed 
within the finite range post form DWBA theory as discussed in section 3.1. The 
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structure inputs were also the same as those given their. The first order results 
have been obtained by using Eq. (81). 

In Fig. (14), we present calculations for the triple differential cross sections 
for the breakup reaction ^^Be + Pb n + ^°Be + Pb, as a function of the 
energy of the ^°Be core (Ec), for four beam energies lying in the range of 5 
MeV/nucleon - 72 MeV/nucleon. To see the postacceleration in a clear way, it 
is very useful to study the cross-section as a function of the core energy. The 
results obtained within the higher order and the first order theories are shown 
by solid and dotted lines, respectively. 

It can be seen from this figure that while for lower beam energies, the higher 
order and first-order results differ considerably from each other, they are almost 
the same for the beam energy of 72 MeV/nucleon. In each case, the first order 
cross sections peak at the energy of the core fragment which corresponds to 
the beam velocity (this value of the core fragment energy will be referred to as 
Ef,„ in the following). In contrast to this, the peaks of the higher order cross 
sections are shifted to energies > E^y for the three lower energies. Only for the 
72 MeV/nucleon beam energy, does the higher order result peak at E^v- This 
shows very clearly that the finite range DWBA model exhibits postacceleration 
for beam energies < 30 MeV/nucleon, while this effect is not present at 72 
MeV/nucleon. Therefore, the higher order effects are minimal for the Coulomb 
breakup of ^^Be at the beam energies > 70 MeV. This result is in agreement 
with those obtained in [1131 1^ . 

In Fig. (15), we compare the results of the first-order and the finite range 
DWBA calculations for the relative energy spectrum of the fragments emitted in 
the breakup reaction of -'^^Be on a ^"^Pb target for the same four beam energies as 
shown in Fig. (14). These cross sections have been obtained by integrating over 
all the allowed values of the angles flc-n- In both the models, the integrations 
over have been carried out between 1° to grazing angle, in the upper 

two figures, and between 5° to grazing angles, in the lower two figures. The 
integrations over angles have been done over all of its kinematically 

allowed values. The dotted and solid lines represent the first-order and the 
higher order results, respectively. 

We notice that while for the beam energy of 72 MeV/nucleon, the higher 
order effects are minimal, they are quite strong for the lower beam energies, 
being largest at the beam energy of 5 MeV/nucleon. This reinforces the point, 
already made in |113l IH2j . that at the beam energy of 72 MeV/nucleon, the 
higher order effects are quite small if both the first order and the higher order 
terms are calculated within the same theory. 

In Fig. (16), we show the same results as in Fig. (15) but for the ^^C in- 
duced reaction on the ^ospj-, i^^rget for the beam energies of 67 MeV/nucleon, 30 
MeV/nucleon and 10 MeV/nucleon. We see that in this case too the higher or- 
der effects are quite weak for the beam energy 67 MeV/nucleon, but appreciable 
for the lower beam energies. 

It may be noted that by comparing the result of the adiabatic model of 
Coulomb breakup reactions which is conceptually different from ours, with that 
of the first order semiclassical perturbation theory of the Coulomb excitation, it 
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Figure 15: The differential cross section as a function of the relative energy of 
the fragments (neutron and ^°Be) emitted in the ^^Be induced breakup reaction 
on a ^"^Pb target at the beam energies of 72 MeV/nucleon, 30 MeV/nucleon, 
10 MeV/nucleon, and 5 MeV/nucleon. The dotted and full lines represent the 
first-order and the finite range DWBA results, respectively. 



has been concluded in 117 that the higher order effects are substantial for these 
reactions even at the beam energies ~ 70 MeV/nucleon. However, one should 
be careful in drawing definite conclusions about the role of the higher order 
effects from such an approach. For a reliable assessment of the contributions of 
the higher order effects, it is essential that both the first order and the higher 
order terms should be calculated within the same theory, as has been done in 

usiEiinsi)- 

Thus, in the post form DWBA theory, the peaks in the triple and double 
differential cross sections vs. core energy spectra, are shifted to energies larger 
than those corresponding to the beam velocity, at the incident energies < 30 
MeV/nucleon. Therefore, postacceleration effects are important at these beam 
energies. On the other hand, at the beam energy ~ 70 MeV/nucleon, the 
corresponding spectra peak at the beam velocity energies, which is consistent 
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Figure 16: The differential cross section as a function of the relative energy of 
the fragments (neutron and ^*C) emitted in the ^^C induced breakup reaction 
on a ^°*Pb target at the beam energies of 67 MeV/nucleon, 30 MeV/nucleon, 
and 10 MeV/nuclcon. The dotted and full lines represent the first-order and 
the finite range DWBA results, respectively. 



with no postacceleration. In contrast to this, the first-order cross sections always 
peak at the beam velocity energy, which is expected as the postacceleration is 
a higher order effect. 

The higher order effects are also found to be quite important in the relative 
energy spectrum of the fragments at beam energies < 30 MeV/nucleon, while 
they are insignificant at the beam energies ~ 70 MeV/nucleon. This suggests 
that the conclusions arrived at in Refs. |90[lll2j . where the data on the relative 
energy spectra of the fragments taken in the breakup of ^^Be and ^^C at beam 
energies ^ 70 MeV, have been analyzed within the first order theory of the 
Coulomb excitation, may not be affected by the higher order effects. 

It should be noted that from an experimental point of view, the postaccel- 
eration effects are not fully clarified (see, e.g., PUI I118L lliyp . Finally, let us 
mention the recent work on the electromagnetic dissociation of unstable neutron- 
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rich oxygen isotopes |12U| . These authors deduce photoneutron cross sections 
from their dissociation measurements. If the neutrons are emitted in a slow 
evaporation process in a later stage of the reaction, the question of postacceler- 
ation is not there. On the other hand, for the light nuclei there is some direct 
neutron emission component and the present kind of theoretical analysis further 
proves the validity of the semiclassical approach used in |120) . 

Postacceleration effects are also of importance for the use of Coulomb disso- 
ciation for the study of radiative capture reactions of astrophysical interest. We 
expect that our present investigations will shed light on questions of postaccel- 
eration and higher order effects in these cases also. 

5 Summary, Conclusions and Future Outlook 

Due to the potential use of the breakup reactions of neutron rich light exotic 
nuclei in extraction the information about the structure of these nuclei, it is es- 
sential to have a full quantal theory of these reactions which involves minimum 
number of parameters. To this end, we have developed a theory of one-neutron 
halo nuclei within the framework of the post form distorted wave Born approx- 
imation. Finite range effects of the core-halo interaction are included in this 
theory which allows the full ground state wave function of the projectile corre- 
sponding to any orbital angular momentum structure to enter into this theory. 
For the one-nucleon halo breakup case, it can treat the Coulomb and nuclear 
breakups as well as their interference terms within a single framework on an 
equal footing. Since this theory uses the post form scattering amplitude, the 
breakup contributions from the entire continuum corresponding to all the multi- 
poles and relative orbital angular momenta of the core-halo system are included 
in it. This theory can account for the postacceleration effects in a unique way. 

Most of the breakup observables are sensitive to the ground state config- 
uration of the projectile. We find that for ^^Be, as- wave configuration 
[-^°Be(0"'") (g) lsi/2v], with a spectroscopic factor of 0.74 for its ground state 
provides best agreement with the experimental data in all the cases. 

We have also performed calculations within an adiabatic model which makes 
the approximation that the strongly excited core-valence particle relative ener- 
gies are small in the Coulomb breakup and also within the approximation of 
Baur and Trautmann which equates the coordinates of the core-target system 
with those of the projectile-target system. Unlike the DWBA, the adiabatic 
model does not use the weak coupling approximation to describe the center of 
mass motion of the fragments with respect to the target. For almost all the 
observables, there is a general agreement between the DWBA and the adiabatic 
model results even in the absolute magnitude. However, the approximation of 
Baur and Trautmann gives results which are very different from those obtained 
within the DWBA and the adiabatic model. 

For the ^^C case, the results for the PMD of ^^C and the relative energy 
spectrum of the n - ^*C system show that the most probable ground state 
configuration of ^^C is [^^C (0"'") (8) lsi/2J^] with a one-neutron separation energy 
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of 530 keV and a spectroscopic factor of 1. Furthermore, the most probable 
configuration for ^^C is a s - wave valence neutron coupled to the ^^C core and 
that for ^^C is a d - wave valence neutron coupled to the ^^C core. Both the 
experimental and the calculated FWHM of the PMD for the ^"'C core in the 
breakup of ^^C are small and they agree well with each other. This provides 
support to the existence of a one- neutron halo structure in ^^C. On the other 
hand, in the case of ^'^C the FWHM of the PMD for the ^^C core is closer 
to that of a stable isotope. Therefore the existence of a halo structure in ^'^C 
appears to be unlikely. Interestingly the one- neutron separation energies of ^^C 
and ^^C are 1.2181 and 0.729 MeV, respectively. So both the binding energy of 
the valence neutron as well as its configuration with respect to the core together 
decide whether a nucleus has halo structure or not. 

For medium mass and heavy target nuclei, the neutron angular distributions 
are dominated by the nuclear and the Coulomb breakup terms at larger and 
smaller angles, respectively. Contributions of the Coulomb-nuclear interference 
terms are non-negligible. They can be as big in magnitude as the pure nuclear 
or the pure Coulomb breakup and have negative or positive sign depending 
upon the angle and energy of the outgoing fragments. For these targets, the 
interference terms help in better description of trends of the experimental data 
even at smaller angles. 

In the parallel momentum distribution of the ^''Be fragment in the breakup 
reaction of ^^Be, the region around the peak of the distribution, which gets sub- 
stantial contributions from forward scattered fragments, is Coulomb dominated, 
while in the wings of the distribution, where contributions come from fragments 
scattered at large angles, the nuclear breakup contributions dominate. Simi- 
larly, the data on relative energy spectra of the fragments (neutron and ^°Be) 
emitted in breakup of ^^Be on a heavy target can not be described properly 
by considering only the pure Coulomb breakup mechanism; inclusion of nuclear 
and Coulomb-nuclear interference terms is necessary. 

In future studies the full quantal theory of the one-neutron halo breakup 
reactions should be apphed to describe the (a, 67) reaction. The data for these 
reactions |121l I122| taken at the Michigan State University are yet to be de- 
scribed within a full quantal reaction model. Since, the relevant partial cross 
sections of the core fragments are essentially inclusive, both elastic and inelastic 
breakup modes will contribute to them. Our theory can be extended to calcu- 
late the latter mode in a straight forward manner. In fact all the ingredients 
required for this extension have already been calculated by us. 

There is also a need to extend the theory to describe the halo breakup at 
higher beam energies for which data have been taken at GSI, Darmstadt. This 
can be achieved by introducing the cikonal expansion of the distorted waves, 
instead of the partial wave expansion as done in section 3.2. 

Our theory should be used to analyze the breakup data of *B on a '^^Ni tar- 
get at the beam energy of 25.8 MeV. At this energy the nuclear breakup effects 
are quite strong and the Coulomb-nuclear interference terms should manifest 
themselves in an important way. It should be noted that in the CDCC model 
analysis of these data, the continuum states corresponding to much larger excita- 
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tion energies and relative orbital angular momenta were required to be included 
in order to get a proper convergence. Our post form breakup theory includes, 
by its very construction, contributions from the entire continuum corresponding 
to all the excitation energies, multipoles and relative orbital angular momenta 
of the core- valence nucleon system. Therefore, a comparison of our calculations 
with these data would be interesting also from the point of view of checking and 
supplementing the corresponding CDCC results. 

The study of the ^Li(a,n)^^B reaction is of crucial importance in determin- 
ing the abundance of ^^B, since nuclides with A > 12 pass through ^^B on 
their way to higher masses. However, the a-capture reaction ^hi{a,n)^^Ti is in 
competition with the *Li(n,7)^Li reaction which turns the reaction flow back 
to lighter elements via ®Li(/3~,i/)^Be(p,a)^Li. Thus the cross sections of these 
reactions are extremely important in predicting the yields of elements with A > 
12. We would like to apply our theory to investigate the breakup reaction ^Li 
— > ^Li + n on a heavy target like ^"^^Pb. It would be interesting to know how 
far the data taken recently at the Michigan State University can be described 
solely by the Coulomb breakup mechanism and what is the role of the nuclear 
breakup effects in these data. This is important as the pure Coulomb breakup 
cross sections for this reaction can be used to extract the cross sections for the 
reverse reaction ®Li(n,7)^Li. 

The authors would like to express their thanks to P. Banerjee, G. Baur, P. 
Danielewicz and H. Lenske for several useful discussions on the present topic. 
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